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Summary 

Analysis and experimental measurement of the 
electromagnetic force loads on the hybrid rotor in a novel 
hybrid magnetic-bearing switched-reluctance motor 
(MBSRM) have been performed. A MBSRM has the 
combined characteristics of a switched-reluctance motor and 
a magnetic bearing. The MBSRM discussed in this report has 
an eight-pole stator and a six-pole hybrid rotor, which is 
composed of circular and scalloped lamination segments. The 
hybrid rotor is levitated using only one set of four stator 
poles, while a second set of four stator poles imparts torque to 
the scalloped portion of the rotor, which is driven in a 
traditional switched reluctance manner by a processor. Static 
torque and radial force analysis were done for rotor poles that 
were oriented to achieve maximum and minimum radial force 
loads on the rotor. The objective is to assess whether simple 
one-dimensional magnetic circuit analysis is sufficient for 
preliminary evaluation of this machine, which may exhibit 
strong three-dimensional electromagnetic field behavior. Two 
magnetic circuit geometries, approximating the complex 
topology of the magnetic fields in and around the hybrid 
rotor, were employed in formulating the electromagnetic 
radial force equations. Reasonable agreement between the 
experimental and the theoretical radial force loads predictions 
was obtained with typical magnetic bearing derating factors 
applied to the predictions. 

Introduction 

There is a need for reliable, fail-safe, robust, compact, low- 
cost electric motors for applications with high temperatures 
or extreme temperature variations; switched-reluctance 
motors possess these characteristics (ref. 1). These motors 
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have been evaluated as high-speed starter-generators (refs. 2 
and 3). However, conventional switched-reluctance motors 
can suffer from undesired vibration due to (1) unbalanced 
lateral forces on the rotor caused by electrical faults, (2) 
mechanical offset of the rotor, and (3) uncontrolled pulsed 
current in one or more shorted coils (ref. 4). A viable solution 
for mitigating mechanical vibration is to suspend the rotor 
magnetically via magnetic levitation. In addition, magnetic 
suspension of the rotor allows the motor to operate at a much 
higher rotational frequency for a prolonged period. This 
benefit is due largely to the elimination of friction, as there is 
no physical contact between the stator and rotor. 
Consequently, motors incorporating magnetic bearings 
perform at higher efficiency than motors incorporating ball 
bearings. Magnetic levitation also obviates the need for a 
lubrication system, which has the added benefit of 
significantly decreasing the weight and complexity of 
turbomachinery mechanisms. 

Methods for simultaneously levitating and rotating a rotor 
within a single stator in switched-reluctance motors have 
been proposed in references 5 and 6. In these motors, the 
technique of using differential stator windings was employed. 
The studies in references 5 and 6 center primarily around the 
use of a main four-pole winding to rotate an eight-pole rotor 
while utilizing a two-pole winding to apply radial force to the 
rotor with all of the 12 stator poles having both windings 
thereon. A variation on this theme was described in 
reference 7, wherein only a single coil on each stator pole (in 
a 12/8 stator-rotor pole configuration) was employed to 
achieve motor-bearing action. 

Self-levitation (also called self-bearing) of motors has been 
achieved for nearly every type of electric motor, but is very 
marginal in performance for switched-reluctance motors with 
low numbers of poles. The motoring technique disclosed in 
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reference 8 will simultaneously levitate and rotate a rotor, not 
only for (18/12) and (12/8) stator-rotor pole combinations, 
but also for (8/6) and (6/4) configurations while employing a 
single set of coils positioned on each stator pole. The 
motoring techniques described in references 5 to 7 are not 
applicable to motors having low stator-rotor pole 
configuration of (8/6) or (6/4) because of the dearth of stator- 
rotor poles in appropriate positions to apply levitating forces. 
The hybrid rotor technique described in reference 8 will 
assure robust bearingless operation in all four aforementioned 
stator-rotor pole configurations. With circular laminations on 
the shaft, levitating force is always assured as the rotor spins. 

A disassembled eight-pole stator, six-pole hybrid rotor 
(motor) is shown in figure 1 . The rotor lamination stack has a 
cylindrical portion of its length (indicated by the left arrow in 
the picture) used for levitation and a scalloped portion 
(indicated by the right arrow) used for both levitation and 
motoring. An algorithm that handles both the levitating and 
motoring functions (ref. 8) controls the currents in the coils 
on the eight-pole stator. During normal operation, two pairs 
of opposing stator poles (at right angles to each other) levitate 
the rotor. The remaining two pairs of stator poles exert torque 
on the six-pole rotor lamination stack to produce rotation. 
The relative lengths of the circular and multipole lamination 
stacks on the rotor can be chosen to tailor the performance of 
the motor for a specific application. For a given overall 
length, increasing the length of the multipole stack relative to 
the circular stack results in an increase in torque relative to 
levitation load capacity and stiffness, and vice versa. The 
motor discussed in this report has a 2.525x1 0” 2 -m- (1-in.-) 
long stator. The rotor has a 2.525xl0” 2 -m- (1-in.-) long 
hybrid segment comprising a 5.883xl0” 3 -m- (0.23 3 -in.-) long 
section of circular laminations for levitation and a 
1.937xl0” 2 -m- (0.767-in.-) long portion having six salient 
poles for rotation. 

This report focuses exclusively on the nonrotating 
magnetic bearing action forces and static torque in a novel 
(8/6) hybrid motor. Three types of magnetic force 
measurements were made in which two diametrically 
opposed rotor poles were aligned with stator poles. In the first 
measurement, only one stator pole was energized, and the 
force on the centered rotor was recorded as a function of 
current. In this situation, the flux through the excited pole 
returns through a number of unexcited poles. In the second 
type of measurement, a proportional-derivative (PD) 
algorithm was used to levitate the nonrotating rotor. The rotor 
was then vertically displaced using a fixture containing a load 
cell, and the stator coil currents and force exerted on the rotor 
were recorded for each displacement value. In the 
corresponding analysis, the measured currents were utilized 
directly, thus making the calculations independent of the 
control law. The third measurement was the determination of 
static torque for rotor angles ranging from alignment to 30° 
from alignment. The question to answer is, “can one 
successfully apply one-dimensional magnetic field analysis to 
the complex three-dimensional hybrid rotor to obtain the 
associated electromagnetic radial forces?” 



Figure 1. — Disassembled eight-pole stator and six-pole 
hybrid rotor of magnetic-bearing switched -reluctance 
motor (MBSRM). 

A list of symbols used in this report is given in the 
appendix to aid the reader. 


Derivation of General Magnetic Force 
Equation 

Figures 2 and 3 depict a portion of the rotor and a portion 
of a stator tooth on which the stator coil is energized to 
produce a magnetic field flux density B. The magnetic field 
interacts with the rotor at a distance g. The magnetic energy 
w is given by the volume integral equation 

*"if fl2dF (i> 


where p 0 is the permeability of free space. Integration over 
only the flux volume V is possible if the magnetic core 
material is operating within its linear range so that the energy 
stored in the core can be neglected compared with that in the 
gap. Referring to figure 2, it can be shown from equation (1) 
that the magnetic radial force on the rotor toward the stator is 


F = - 


d w 

dg 


B z 

2 l^o 


11 iw = ~EA 


( 2 ) 


where dg is the differential gap length, L c is the stack length 
of the circular laminations, is the width of a stator tooth, E 
is the general energy density between a stator pole and the 
rotor, and A is the general common cross-sectional area 
between a stator pole and the rotor surface. Similarly, the 
instantaneous torque on the rotor (fig. 3) at the scalloped 
segment is given by 


T = - 


d w 
dO 


B 2 T S 

L gr partial overlap 

< 2 M 0 

0 aligned, no overlap 


( 3 ) 
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Figure 2,— Side view of circular rotor and stator tooth 
segment of 8/6 MBSRM, where B is magnetic flux 
density, g is gap length, V is tooth width, ris outer 
radius of rotor laminations, and V is flux volume. 



Figure 3.~Side view sketch of scalloped rotor-stator tooth seg- 
ment, where g is gap length, / is rotor tooth width, r is outer 
radius of rotor laminations, and dV is differential flux volume. 


where d0 is the differential angular displacement, L s is the 
stack length of the scalloped laminations, g is the gap length, 
and r is the outer radius of scalloped lamination section of the 
rotor. The mathematical form and value of the radial B field 
impinging on the nonrotating asymmetric rotor will be 
determined next. 

Nonlevitated Magnetic Equivalent 
Circuit Modeling for Centered Rotor 

Figure 4 is an axial view of the motor. 

Circular Lamination Force Contribution 

A complex three-dimensional magnetic field topology 
exists in and around the rotor because of its hybrid design. 
However, the analysis will be restricted to a one-dimensional 
approach in order to reduce the mathematical complexity of 
the static radial force derivation. Each rotor segment (circular 
lamination stack and scalloped lamination stack) will be 
treated as an independent entity, and accordingly, we develop 
a one-dimensional magnetic radial force description for each 
rotor segment. The magnetic radial forces will then be 
summed to obtain the total static magnetic radial force on the 
rotor. 

Figure 5 depicts the equivalent two-dimensional flux 
distribution in the rotor’s circular laminations and the eight 
stator poles. Each stator pole is at a distance g 0 (gap length) 
from the circular lamination stack and each pole shares a 
common cross-sectional area A c with the circular laminations. 
Figure 6 is the corresponding magnetic equivalent circuit of 
figure 5. Both figures will be employed in deriving the static 
magnetic radial force equations. 


Pole 5 



pole 1 

Figure 4. — Axial view of the 8/6 hybrid MBSRM with rotor 
and stator pole 1 aligned. 

The common cross-sectional area A c between the circular 
lamination stack and each stator pole can be expressed as 

A£=A C , k - 1...8 (4) 

where A c - L c l tw . The flux at stator pole 1 (|>f , produced by 
energizing coil 1, is distributed equally among the other 
seven poles. The flux distribution obeys the relation 

k = 2 
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Figure 5, — Equivalent two-dimensional rotor-stator flux 
distribution for nonlevitated 8/6 hybrid MBSRM, where 
for the circular lamination ^ is magnetic flux between 
stator pole k and rotor lamination stack, A £ is 
cross-sectional area between stator pole k and rotor 
lamination stack, and go is gap length. 


and noting that mmf= ni , where n is the number of turns in 
the coil and i is the current in the coil, we can write the flux 
expression 



( 10 ) 


Substituting equation (8) into (10) yields 



\i 0 A c ni 

So 


( 11 ) 


from which we obtain the flux density at pole 1 : 
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A c go , 


( 12 ) 
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Figure 6. — Magnetic equivalent circuit for nonlevitated 
8/6 hybrid MBSRM, where for circular lamination 
is reluctance and is magnetic flux between stator 
pole k and rotor lamination stack. 


and for a rotor that is radially centered, the reluctances 
between the stator poles k and the rotor are related by 


1| 

49 


( 2 2-2 \ 

p 0 n 1 

2l 

,64 ) 

UoJ 

K So ) 


Substituting equation (12) into an expression similar to 
equation (2), we can write the magnetic radial force as 


(13) 


where the B field is considered uniform. Equation (13) is the 
static magnetic radial force due to pole 1 (fig. 4). The 
oppositely directed radial forces produced by poles 3 and 7 
will sum to zero. Likewise, the resultant magnetic force 
produced by poles 2, 4, 6, and 8 is zero. Hence, only the 
radial force at pole 5 will contribute to the resultant force on 
the rotor. Because of the flux symmetry in the rotor and 
stator, we can write the magnetic field at pole 5 as 


K = K+ 1. * = i-7 


( 6 ) 


b c 5 =-b? 

5 7 


(14) 


The total reluctance for the flux passing through the powered 
pole is 


n c — i) c i n c 
lYj’ — Ivy + 4V|| 


(7) 


and substituting equation (12) into (14) yields 


B | =- 

J o 


So ) 


(15) 


where is in series with the parallel reluctances 

R \ , , Rl , Rj , and R$ . Further, it can be shown 

that equation (7) has the form 


R 


c 

T 


8 

7 


f \ 

So 

j 


( 8 ) 


From the basic magnetic circuit equation, the magnetomotive 
force 


The magnetic radial force at pole 5 can now be written as 


If 

'n 


( 2 2 -2\ 

p 0 n 1 

21 

,64 J 

UoJ 

So J 


(16) 


and the net downward magnetic force produced by poles 1 
and 5 on the circular laminations is given as 


mmf = (])/? 


(9) 


Ket=F l 


(17) 
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Substituting equations (13) and (16) into (17) yields the net 
downward magnetic force on the circular lamination stack; 


F = 

net 


3 f A c yi 0 n 2 ^ 


S o 


(18) 


Substituting the actual value of the constants gives 

F c net = 5.037 r (19) 


We now determine the static magnetic radial force on the 
scalloped laminations by employing some aspects of the 
previous discussion. 



Figure 7. — Equivalent two-dimensional rotor-stator 
flux distribution for nonlevitated 8/6 hybrid MBSRM, 
where for scalloped lamination g 0 is gap length, ^ 
is magnetic flux between stator pole k and rotor 
lamination stack, and A s k is cross-sectional area 
between stator pole k and rotor lamination stack. 


Scalloped Lamination Force Contribution 

Figure 7 depicts the equivalent two-dimensional flux 
distribution in the rotor scalloped laminations and their 
adjacent stator poles. Each rotor pole is at a nominal distance 
go from its adjacent stator poles. Four stator poles — 2, 4, 6, 
and 8 (figs. 4 and 7) — share a common cross-sectional area 

■jA s with four adjacent rotor poles. Figure 8 is the 

corresponding magnetic equivalent circuit of figure 7. 
Additionally, in the spirit of rough approximation, we neglect 
any fringing effects from the rotor to stator poles 3 and 7. The 
expected large reluctance between the stator poles and the 
scalloped cavities will result in a relatively small return flux. 
However, in later sections analyzing the unaligned centered 
rotor and the aligned levitated rotor, the fringing effect will 
be included to give a more complete description of the 
levitation force on the rotor. 

Observation of figure 4 reveals that only the stator poles 
labeled as 1,2, 4, 5, 6, and 8 provide significant paths for the 
circulating magnetic flux, and each pole 2, 4, 6, and 8 has an 
equivalent, effective common stator-rotor area that is one- 
third that of stator pole 1 or 5. The common cross-sectional 
area between the rotor scalloped lamination stack and a stator 
pole can be expressed as 


4 = A s 5 = A s 


( 20 ) 


w < 

c < 

c < 

: R| : 

; sR| : 

> r| < 
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: «a < 
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, , 
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Figure 8. — Magnetic equivalent circuit for 
nonlevitated 8/6 hybrid MBSRM, where 
for scalloped lamination is reluctance 
and is magnetic flux between stator 
pole k and rotor lamination stack. 

The total reluctance in the magnetic circuit can be written as 

R s t = Rf + Rfi (24) 

where R( is in series with the parallel reluctances 

R 2 , R 4 , ^5 , R$ , and R$ . The resultant parallel reluctance is 

thus obtained by 


*11 = 


1 



(25) 


and 


4k=\4=\A s , k = 1...4 (21) 

The reluctance between the rotor and stator poles for the 
radially centered rotor is hence given as 


Substituting equation (25) into (24) and using equations (22) 
and (23), it can be shown that the total reluctance at pole 1 is 


R 


s 

T 


10 

7 


f \ 

go 

y 


(26) 


and 


go J^s 

Po4 s 


go _ 3 gp _ 

u 4 Mf 1 

^ 3 


Ri=R s 5 and R'i k = /?f , k = 2..A 


Using the magnetomotive force argument, we can derive 
^ ' equation (27) 



f A s A 

p 0 A x m 

v g ° J 


(27) 


(23) The flux density extant at pole 1 and its adjacent salient rotor 
pole is thus given as 
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1_ 

10 


f \x {] ni ' 

go 


(28) 


Substituting equation (28) into an expression similar to 
equation (2), the magnetic force at pole 1 is given as 


Fi 


f— 1 

U) 

( 2 2-2 ) 

p 0 n l 

UooJ 

UoJ 

go J 


(29) 


77 — J7 C _i_ J7 S 

r t{net ) r net " r r net 


Substituting equations (19) and (36) into (37) yields 


Ff(net) = !3- 88 / 


(38) 


The fraction of force contributed by the circular lamination 
segment is 


Equation (29) represents the magnetic force, due to stator 
pole 1 , on the aligned salient pole of the rotor. The resultant 
magnetic force on the rotor, produced by poles 2 and 8, 
cancels the resultant magnetic force produced by poles 4 and 
6. Hence, only the magnetic force at pole 5 will contribute to 
the resultant force on the rotor. The flux at pole 5 can be 
expressed as 


^5 “ ( 30 ) 

k = 1 

The flux relation at the poles can be written as 


4>2*=<l>2, * = 2...4 (31) 


and making the appropriate substitutions it can be shown that 



f AS A 

p 0 A 5 m 

V SO , 


(32) 


Using equation (32), the magnetic force between pole 5 and 
its adjacent salient rotor pole is 


w 

— 1 

U) 

Vo 

21 

,100 J 

UoJ 

go J 


(33) 


The net downward magnetic force produced by poles 1 and 5 
on their corresponding salient rotor poles is 


F s 

r net 


■ F{ -F 5 s 


(34) 


Substituting equations (29) and (33) into (34) yields the force 
expression 


F s - — 

± net ^ 


f js 2 ^ 

A p 0 n 

si 


(35) 


where Equation (35) represents the net 

downward magnetic force on the scalloped lamination stack. 
Substituting the actual values of the constants gives 

F^= 8.845 i 2 (36) 

The total net downward magnetic force on the centered 
nonlevitated hybrid rotor is the sum of the net downward 
magnetic force on the circular lamination stack and the net 
downward magnetic force on the scalloped lamination stack: 


5.037 i --o.36 = 36percen t (39) 

F t{ net) 13.88 i 2 

and the fraction by the scalloped segment is 

^ net - z ~ 0.64 = 64 percent (40) 

F,(net) 13.88 i 2 

The circular lamination stack experiences approximately 
36 percent of the total net downward magnetic force, while 
the scalloped lamination experience approximately 64 percent 
of said force. The relative stack lengths of the circular and 
scalloped laminations are L c = 0.23 and L s = 0.77, 
respectively. Note that the forces in equations (39) and (40) 
are not proportional to the stack lengths. This discrepancy is 
due to the differences in return paths of the flux in the hybrid 
rotor. The following section outlines the experimental 
procedure for obtaining the nonlevitated magnetic radial 
force on the rotor. 

Experimental Centered-Force Procedure 

The experimental nonlevitated magnetic radial force on the 
rotor was obtained by mounting the rotor between two load 
cells (see fig. 9) while ensuring that the rotor’s axis was 
radially centered in the stator. A salient rotor pole was 
aligned with stator pole 1 (see fig. 4). The preload voltages of 
the cells were then zeroed to eliminate the preload weight 
from the data. Stator pole 1 was energized, and the resulting 
magnetic force on the rotor produced a larger reading on the 
lower cell and a smaller reading on the upper cell. The 
difference in the readings gave the net force on the rotor. This 
procedure was repeated for several values of coil current 
ranging from 0.982 to 2.027 A. The procedure was repeated 
with pole 5 as the sole energized pole. 



Figure 9,— Experimental setup used in obtaining 
nonlevitation magnetic radial forces for 8/6 
hybrid MBSRM. 
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Experimental and Theoretical Centered-Force 
Comparison 

The experimental and theoretical magnetic radial force 
data plot are presented in figure 10. Because of flux leakage 
and field fringing effects at the stator and rotor pole edges, 
uncertainties in geometry, rotor-stator alignment, and 
material properties, we must apply a derating factor to 
equation (38). Reference 9 presents a detailed treatment of 
the derating factor, wherein a factor of 0.88 was determined 
for one particular active magnetic bearing. In contrast, a 
graph in reference 10 shows a derating factor of 0.6 for our 
gap go (5xl0 -4 m). For this hybrid rotor configuration, which 
clearly is different from that of a pure magnetic bearing, a 
derating of 0.87 gives the closest agreement with the 
experimental results (see eq. (41) and fig. 10): 

F t{net) =(13.88x0.87) i 2 (41) 

Figure 10 depicts the experimental least-squares-fit 
magnetic force curves for the upper and lower poles, and the 
theoretical magnetic force curve. In the following section, the 
nonlevitated unaligned rotor-stator radial force equations are 
derived for the hybrid rotor. 

Nonlevitated Unaligned Rotor Analysis 

The analysis presented in previous sections clearly 
demonstrates that one can make a reasonably accurate 
prediction of electromagnet forces on a radially centered 
hybrid rotor using two-dimensional visualization and one- 
dimensional electromagnetic force analysis. With this fact 
established, we will estimate, without experimental 
verification, the radial force on the unaligned radially 
centered rotor, that is, a rotor that is rotated so that the 
motoring section has an unfavorable misalignment with the 
powered stator pole 1, as depicted in figure 11. Using a force 
symmetry argument similar to that presented in earlier 
sections, it can be shown that only poles 1 and 5 are relevant 
in determining the net force on the hybrid rotor. The analysis 
begins by approximating the reluctance between pole 1 or 5 
and the adjacent scallop cavities: 

Rf = = R s 5 (42) 


where g a = 


gb+Sc 


is the average stator-rotor gap length 


at stator pole 1 or 5 (fig. 11), g b is the distance between a 
comer edge of the salient stator pole 1 or 5 and either of its 
adjacent salient rotor pole corner edges, and g c is the gap 
length between the midpoint of either pole face and the base 


of its adjacent scallop cavity. Here, A a = 


' 4 + 4 ' 


is the 


average of the stator pole face area A b and its adjacent 
scallop cavity area A s c . The face area of pole 1 is 



Figure 10. — Nonlevitation magnetic radial force curves for 
8/6 hybrid MBSRM. Experimental least-squares-fit curves 
(dotted lines) for upper and lower poles and theoretical 
curve (solid line). 
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Figure 11 . — Axial view of 8/6 hybrid MBSRM with rotor and 
stator pole 1 unaligned. 

A]!, - A s - 2x 1CT 4 m 2 , and the surface area inside a scallop 

cavity was determined to be A s c - 4.35x1 CT 4 m 2 . 
Substituting these values into the average stator-rotor gap 
length and area expressions result in g a = 6.25x1 0” 3 m and A a 

= 3.17xl0” 4 m 2 . We can rewrite A a in terms of A s : 

A a - l-59A s (43) 


and g a in terms of g 0 : 


ga =12.5g 0 (44) 

It can be shown by substituting equations (43) and (44) into 
(42) that 


Rf =R s 5 = 7MR s 0 = 1.57xl0 7 H -1 (45) 
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where 4 is the nominal reluctance at the aligned scalloped 
segment. Here, the estimated reluctance at the unaligned pole 
is almost eight times that of an aligned pole. Note that this 
value of the reluctance applies to both poles 1 and 5. The 
remaining pole reluctances are 

R S _gQ_, k = 2. ..4, 6.. .8 


We know that 


F s =- — 


= 0.35/" 


The next task is to determine the flux fraction at pole 5, 
which will allow computation of the magnetic force at pole 5. 
The difference between this force and the force at pole 1 will 
give the net magnetic force on the unaligned hybrid rotor. 

The fraction of the total flux passing through pole 5 is 
determined by its reluctance ratio, the combined parallel 
reluctance of all the return poles divided by the reluctance of 
pole 5: 


4 =4 =A S 


r> V N _ 

^5 ratio „ 


4 0.29 RZ 


■ - 0.037 


4 7.864 


n s — n s _ n s 
K 3 ~ K 7 ~ K 0 


The flux at pole 5 is now obtained by employing 
equation (53): 


jjs So ~*So ^2 jjs 

K 2 ~ 7 7 \ “ “ “ ^ K 0 


A s \i 0 A> 


where 4fr = 4 , k = 2... 4 . The total parallel reluctance due 
to poles 2, 3, 4, 5, 6, 7, and 8 is given as 


or simply 


*5= *5 ratio¥\ = 0.037 X — 
Rt, 


<\> s 5 =4.86x1 0~ 7 i 


The magnetic radial force at pole 5 is 


4 2 1 

1 1 

uS uS uS 

K 2 K 3 K 5 


and substituting equations (45), (48), and (49) into (50) gives 


- = 4.71xlO~V 


The net downward magnetic radial force on the unaligned 
scalloped section of the rotor due to poles 1 and 5 is 


4 = 0.294 


77 s — J7 S 77 s 

t u net ~ ~ ^5 


Hence, the total reluctance of the magnetic circuit through 
pole 1 is 

= R[ + 4 = 7.864 + 0.294 = 8.154 (52) 

The reluctance value derived in equation (52) is used later in 
estimating the levitation load on the rotor. Substituting 

4 = = 1 .99 x 10 6 H -1 into equation (52) gives 

4i = 1.63x1 0 7 H -1 . The total flux at pole 1 has the form 


Substituting equations (55) and (59) into (60) produces 


K net ~ 0-35/" 


The ratio of the net downward magnetic force on the 
unaligned scalloped stack to that on the aligned scalloped 


stack F net is 


■ - 0.039 = 3.9percent 


(|)'i =-^- = (6.15xl0“ 8 H)w 


from which we can determine the flux density: 


B s x = = 6-15x1 ^ — ni = (6.63 x 10“ 2 H/m 2 )i (54) 

1 A s 2xl0 -4 m 2 J ' > 


Using equation (2), the radial force at pole 1 is thus 


Note that the unaligned scallop magnetic force load is only 
3.9 percent of the magnetic force load produced by the 
aligned scallop. This is a reasonable result, considering the 
large reluctances of the scallop cavities located at poles 1 and 
5. Taking into account the contributions of the circular 
laminations, the total net downward magnetic radial force on 
the unaligned rotor is 


77 — 77 c _i_ 77 s 

r t(u _ net ) r net ^ r u _ net 


NAS A/TP— 2008-2 14818 


8 



We know from equations (19) and (36) that the aligned net 
radial force on the circular laminations and scallops are 

Fn et ~ 5.037/ 2 and F„ et - 8.845 i 2 , respectively, from which 
the total net downward magnetic force on the aligned rotor is 


P — /7 C 4- 

r t(net ) r net " r r net 


The ratio of the total net downwards magnetic force of the 
unaligned rotor to the total net downwards magnetic force on 
the aligned rotor is 


F t(u_na) _ (5.037 + 0.348)/ 2 
F t(net) (5.037 + 8.845)/ 2 


- 0.39 = 39 percent (65) 


Equations (62) and (65) highlight the importance of the 
circular lamination stack in this hybrid motor; that is, the load 
capacity in the unaligned orientation is substantially 
enhanced by the circular lamination segment (this is critical 
for maintaining rotor stability during the motoring operation). 
In subsequent sections, the general levitation force equations 
are derived for the hybrid rotor that is displaced vertically 
from its radially centered position within the stator. 

Levitation Magnetic Equivalent Circuit 
Modeling 

The analysis begins by employing the magnetic equivalent 
circuit depicted in figure 12, which relates to the magnetically 
levitated state of the nonrotating, aligned, vertically displaced 
hybrid rotor. The rotor was magnetically levitated using 
stator poles 1, 3, 5, and 7 (fig. 4). 

Levitation Flux Derivation 

We designate the coil currents as i u z 3 , z 5 , and z 7 for the 
respective stator poles 1, 3, 5, and 7. Each coil has n turns. 
These currents generate the magnetic fluxes 4>i, 4> 3 , 4> 5 , 4> 7 , § 2 ,s, 
and (j>4 ?6 , where § 2 ,s and (]>4,6 represent the return fluxes due to 
pole pairs (2 and 8) and (4 and 6) on either side of poles 1 
and 5, respectively. Because multiple poles were energized 
during levitation, we will use the permeance expression 
instead of the reluctance term in order to facilitate easy 
mathematical manipulation; hence, by referencing figure 12, 
equations (66a) to (66e) are written by applying the magnetic 
circuit equation to various closed paths (all of which go 
through pole 1, which has a permeance Pi): 

(b 3 (bi 

— - nU - — - + nu (66a) 


r Z 1 

r'3 

r z 5 

r*7 



<j 




> Pi : 

r < 

P3 : 
< 


| pj < 

5 < 

P 2 ,%i 
> < 

A 1 ! 

, <t>3 . 

, 4*5 

.. ^ ! 

, *1*2.8 , 


Figure 12 — Magnetic equivalent circuit for 
levitated 8/6 hybrid MBSRM, where is 
current in coil k, P ^ is permanence at 
stator pole k , and ^ is magnetic flux at 
stator pole k. 

— = — + ^i (66d) 

^2,8 M 

4*4, 6 4*1 , . x 

—— = — + nh (66e) 

M,6 P 

Conservation of the magnetic flux assures that the sum of the 
magnetic fluxes is equal to zero, hence 

4>1 + 4*3 + 4*5 + 4*7 + 4*2,8 + 4*4,6 = ^ (66f) 

The solution of equations (66a) to (66f) is given by 


P r 

( l>i = 


+(■^3 +P 5 + P 7 + P 2> 8 + ^4 6 )(| ] / 


‘ta - F[ p \h +p 5H - Pi h 


+ (P\ + P 5 + P] + ^2,8 + A,6 )'3 ] 1 


^5 - 77" [^'3 - P\i\ + P]ij 


+ (P\ + P}+Pj+ ^ 2,8 + P 4 , 6 ) ; 5 ] ' 


( l*7 - ~rr[ p \h ~ + P 5h 


+(4 + Pi + ^5 + ^2,8 + A, 6 )h J n 


^2,8 ~~ [-^ 1*1 - Pjh + P 5 I 5 ~ Pjij] n (67e) 


4*5 -4*i 

— + HU = — +ZZZ1 


4*4,6 “ - ^3 Z 3 + ^5 Z 5 “ ^7 Z 7 ] n (67f) 


4*7 - 4^1 

nin = — L + nu 

n ' n 1 


where for notational convenience we define the sum of 
permeances as 

P T = /} + P 3 + P 5 + P 7 + P 2 8 + P 4 6 (67g) 


NAS A/TP— 2008-2 14818 


9 



by employing equation (2), hence 


Equations (67a), (67c), (67e), (67f), and (67g) are used to 
derive the electromagnetic levitation force equations for the 
rotor displaced vertically inside the stator. In subsequent 
sections, separate magnetic force equations for the circular 
and scalloped lamination segments are derived. 

Levitation Force on the Circular Laminations 

The common cross-sectional area between the circular 
laminations and pole 1 or 5 is given as A c = L c l tw . The 
rotor’s vertical displacement x from its nominal position 
affects the gap lengths g k , hence 


g\ =£o +x 

at stator pole 1, and 

(68) 

£5 =£o-* 

(69) 

at stator pole 5. In addition, the gaps at the return flux poles 2 
and 8 can be written as 

£ 2,8 =go+* cos ( 45 °) 

(70) 

and for poles 4 and 6 


£4,6 =go-* cos ( 45 °) 

(71) 

where 45° is the angle between adjacent stator poles. The 
approximate permeances between the poles and the circular 
lamination stack are 

If = pc = pc = M‘ 
£o 

(72) 

pc _ 11 c. pC _ 

r k - , k - 1 , 3 , r 2,8 - 

Sk ^ 2,8 

andP 4 c 6 = 2 ^° AC . 

g 4,6 

i C 

? 

(73) 

For the circular laminations, equations (67c) and (67g) are 
utilized to obtain the flux at pole 5, thus 

- ~~ \_ p 3 h ~ P\h + p i h 
'r 

(74) 

+ (^pc +P c +P c + pc s+ pc 6 y i 5 

n 

The equations for the magnetic field and energy density at 
pole 5 are, respectively, 

4 

(75a) 

EH*)- 1 * 1 ™ 2 

2|4q 

(75b) 


The magnetic force on the rotor due to pole 5 is then obtained 


F 5 c =-E c 5 (x)A c (76) 

where B is considered uniform. Using equation (67f), the 
return magnetic flux equation for poles 4 and 6 is 



~ F ih + F £ z 5 ~ F i h^) n 


(77) 


where P£ is the sum of permeances at the circular 

laminations. The equations for the magnetic field and energy 
density at poles 4 and 6 are, respectively, 


B C 4, 6 (X) = 


4*4,6 O) 
A c 


(78a) 


E C a,6 (*) = 


[*4,6 «] 2 

2^o 


and the magnetic force due to poles 4 and 6 is 


(78b) 


* 4 , 6 =- 


-4,6 


(x)2 A c 


(79) 


The total magnetic force due to pole 5 and its adjacent poles 4 
and 6 is thus 


F T (upper) ~ F 5 + * 4,6 ( 80 ) 

Similarly, we can calculate the net downward magnetic force 
created by poles 1, 2, and 8. The magnetic force at pole 1 is 

Ff=-EZ{x)A c (81) 

and the magnetic force due to the return flux poles 2 and 8 is 

Fl % =-E c ^{x)2A c (82) 

The total magnetic force due to pole 1 and its adjacent poles 2 
and 8 is thus 


F T (lower) ~ F \ + * 2,8 ( 83 ) 

The net magnetic levitation force on the circular laminations, 
due to poles 1, 2, 4, 5, 6, and 8 can now be written: 

F T(net_lev) ~~ F T (upper) F T (lower) ( 84 ") 

Note that for the vertically displaced rotor, the net force on 
the circular lamination stack, due to poles 3 and 7, is 
identically zero. 

Levitation Force on the Scalloped Laminations 

We now determine the net force on the scalloped 
laminations. The common cross-sectional area between the 
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scalloped laminations and pole 1 or 5 is given as A s - 
where V is the scalloped rotor lamination stack length. For 
the return poles 2 or 8 and 4 or 6, the common cross-sectional 

areas are ^(2,s) = F l tw(2,s) an( * A ( 4 , 6 ) = F l tw( 4,6) 

respectively, where the corresponding tooth widths are 

_2 

W,8) = — p xsin(51.6°) (85) 

and 

_2 

L(4,6) = LQ3 X 3 10 + x sin(5 1 .6°) (86) 

Here, 51.6° is the angle between the midpoint of stator pole 1 
and either of the midpoints of the overlap of its adjacent poles 
2 or 8 with their adjacent rotor poles, and similarly for stator 
pole 5 and its adjacent poles 4 or 6. The approximate fringing 

permeance P$ or Pj at pole 3 or 7, respectively, is similar in 
value to that used in equation (53), hence 

P{ =p 7 s = J_ = 6.1509x10“ 8 H (87) 

r t\ 

The remaining permeances are 


4*4,6 ~ _ h + F s z 5 _ F i h ) n (92) 

P T 

The equations for the magnetic field and energy density at 
poles 4 and 6 are, respectively, 


R s ^4,6 

4 ' 6 A s 

a 4,6 

(93a) 

c s , , [«4,6(X)] 2 

^4,6 W= 0 

2o 0 

(93b) 


and the equation for the magnetic force due to poles 4 and 6 
is now given as 

Fl 6 =-E{ 6 (x)2A s (94) 

where the gap expression for differentiation is 
g 4 6 = g 0 _xcos (51-6°) . 51.6° is the angle between the 

midpoint of stator pole 5 and either of the midpoints of the 
overlap of its adjacent poles 4 or 6 with their adjacent rotor 
pole. The equation for the total magnetic force due to pole 5 
and its adjacent poles 4 and 6 on the scalloped laminations is 
thus 


% 


Sk 


k = 1,5; 


P s — 

r 2,8 ~ 


2pp^2,8 . 
& 2,8 


and P 4 (y — 


2Fo^4,6 

&4,6 


( 88 ) 


For the scalloped laminations, equations (67c) and (67g) are 
utilized to obtain the flux at pole 5, thus 


F T (upper) ~ F 5 + ^4,6 ( 95 ) 

Similarly, we can calculate the net magnetic force on the 
scallop laminations created by poles 1, 2, and 8. The 
magnetic force at pole 1 is 


F{ = -Ei (x)A s (96) 


and the magnetic force due to the return flux poles 2 and 8 is 



F 3 h ~ F \h + F i z 7 


+ + P^ + Pi + P{ 8 



n 


(89) 


F is=~ F 2,s(x)2A s (97) 

The total magnetic force due to pole 1 and its adjacent poles 2 
and 8 is hence 


The equations for the magnetic field and energy density, at 
pole 5 are, respectively, 



E s 5 (x) = 


2|x 0 


(90a) 

(90b) 


The magnetic force on the scalloped laminations, due to pole 
5, is thus given as 


F T (lower) ~ F \ + F 2,S ( 98 ) 

The net magnetic levitation force on the scalloped 
laminations due to poles 1, 2, 4, 5, 6, and 8 is thus 

F T(net_lev) ~ F T (upper) ^ F T (lower) (99) 

Note that for the vertically displaced rotor, the net force on 
the scalloped lamination stack due to poles 3 and 7 is 
identically zero. The total net magnetic levitating force on the 
vertically displaced hybrid rotor can now be written: 


F 5 s =~E s 5 (x)A s (91) 


r T(net _lev) 


J7 C 

r T(net _lev) 


+ F, 


T (net _lev) 


( 100 ) 


Using equation (67f), the return magnetic flux equation for 
poles 4 and 6 is 


Clearly, the motor becomes a pure magnetic bearing if solid 
rings replace the scallops. In this instance equation (100) will 
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have only the F^ net lev ^ term, in which A c is replaced by a 

value that is representative of total face area of a stator pole. 
The experimental measurements and their comparisons with 
the predictions of equation (100) are presented in the 
following sections. 

Experimental Levitation Force Procedure 

The experimental setup is similar to that displayed in 
figure 9 for the nonlevitated configuration, except that the 
rotor is constrained to move in the vertical plane only. The 
hybrid rotor was configured as shown in figure 4 and then 
magnetically levitated, without rotation, using stator poles 1 , 
3, 5, and 7 and a standard PD control law (see ref. 8). A 
vernier device supporting a load cell, which in turn is 
attached to the rotor via a wire, effects vertical motion. The 
rotor was then preloaded to ensure tension in the load cell for 
both positive and negative vertical displacements and the load 
cell voltage zeroed, thereby eliminating the combined preload 
and rotor weight from the data. Fifteen position data points, 
along with their four coil currents q, z 2 , z 3 , and z 4 , were 
recorded. 

Experimental and Theoretical Levitation Force 
Comparison 

The experimental rotor displacements discussed above 
along with their associated coil currents and the levitation 
force equations derived in (100) were programmed into 
Mathematica (Wolfram Research), which generated the 
predicted force data points. 

Because of (1) flux leakage and field fringing effects at the 
stator and rotor pole edges, (2) experimental errors related to 
the forced rotor displacement, and (3) uncertainties in rotor- 
stator geometry, alignment, and material properties, a 
derating must be applied to the predicted data set. A derating 
of 0.72 gave the closest agreement to the experimental result. 
The absolute values of the experimental and derated 
theoretical data points were least squares fitted using Excel 
(Microsoft Corporation); both are presented in figure 13. 

The static torque analysis for the hybrid rotor is discussed 
next. 


Torque Analysis 

An algorithm that handles both the levitating and motoring 
functions controls the currents in the stator coils. Two pairs 
of opposing stator poles (at right angles to each other) levitate 
the rotor. The remaining two pairs of stator poles exert torque 
on the six-pole rotor lamination stack to produce rotation. 
Figure 14 is a plot of the experimental static torque obtained 
with the aid of a Magtrol dynamometer. Using equation (3), 
the predicted instantaneous torque, generated by two opposite 
stator poles during partial overlap, is given as 


T = 


\i 0 n 2 L s r 

go 


( 101 ) 



Figure 13. — Least-squares-fit experimental (dotted line) 
and theoretical (solid line) levitation magnetic radial 
force curves for 8/6 hybrid MBSRM. 



Figure 14. —Experimental torque graph for 8/6 hybrid MBSRM 
showing torque variation on rotor poles from —4° to 35.5°. 


For a current of 2 A, the predicted instantaneous torque is 
0.227 Nm. The average predicted torque is obtained by taking 
the product of the rotor theoretical duty cycle (22.5°/60°) and 
0.227, which is 0.085 Nm. The average experimental torque, 
produced by the same 2 A current, is similarly obtained to be 
0.058 Nm, using the experimental duty cycle (30760°). It is 
observed that the average experimental torque is 
approximately 0.68 times the average predicted torque. We 
surmise that this lower experimental average torque value 
may be due to significant fringing losses at the poles and the 
influence of the nearby circular laminations. 


Concluding Remarks 

The nonrotating electromagnetic radial forces within a 
novel eight- stator-pole, six-rotor-pole magnetic-bearing 
switched-reluctance motor were examined theoretically and 
experimentally. The circular lamination stack and the 
scalloped lamination stack were treated as independent 
entities, and one-dimensional magnetic force equations were 
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developed for each. The total net magnetic forces on the rotor 
were obtained by summing the net magnetic forces on each 
segment and then applying derating factors of 0.87 for the 
nonlevitated case and 0.72 for the levitated case. These 
derating values are within the range of those discussed in 
previous studies. 

The variation in the radial force derating factors is due, in 
part, to the disparate mode of operation; that is, forced 
displacement of the levitated rotor versus nonlevitation. 
Forced displacement of the levitated rotor introduces an 
additional source of error. A derating factor of 0.68 was 
determined for the average predicted torque. 

In addition, it was demonstrated analytically that the load 
capacity of the hybrid rotor in the unaligned orientation is 
substantially enhanced by the presence of a circular 
lamination segment (this is critical for maintaining rotor 
stability during the motoring operation). Based on the 
analysis, we conclude that two-dimensional visualization and 
one-dimensional magnetic field analysis can be used for 
estimating, preliminarily, the electromagnetic levitation 
forces on the hybrid rotor with a derating factor of about 0.8 
applied to the calculated force. 
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Appendix 

A General common cross-sectional area between a 

stator pole and rotor surface (m 2 ) 

A c Common cross-sectional area between a stator 

pole and circular lamination stack (m 2 ) 

A a Average of A b and A s c (m 2 ) 

A£ Face area of stator pole 1 defined by length of 

scalloped segment, 2x1 0“ 4 m 2 

A s c Surface area inside scallop cavity, rO L\ 

4.35X10 -4 m 2 

A c k Common cross-sectional area between stator 
pole k and circular lamination stack (m 2 ) 

A s k Common cross-sectional area between stator 
pole k and scalloped lamination stack (m 2 ) 

A s Common cross-sectional area between aligned 

rotor and stator pole (m 2 ) 

B General flux density (T) 

B c k Flux density between stator pole k and circular 
rotor lamination stack (T) 

B s k Flux density between stator pole k and scalloped 

rotor lamination stack (T) 

E General energy density between a stator pole 

and rotor surface (J/m 2 ) 

El Energy density between stator pole k and 

circular rotor laminations (J/m 2 ) 

E s k Energy density between stator pole k and 

scalloped rotor laminations (J/m 2 ) 

F k Magnetic force between circular lamination 

stack and stator pole k (N) 

F k Magnetic force between scalloped lamination 

stack and stator pole k (N) 

F Magnetic force (N) 

F„ et Net downward magnetic force on circular rotor 

lamination stack (N) 

F* et Net downward magnetic force on scalloped 

rotor lamination stack (N) 

Ft (lower) Magnetic force on circular rotor lamination due 
to stator poles 1, 2, and 8 (N) 

Ft (lower) Magnetic force on scalloped rotor lamination 

due to stator poles 1,2, and 8 (N) 


-Symbols 

E T(netjev) Total net magnetic levitating force on the 

rotor (N) 

Ft (net lev) Total net ma g ne tic levitation force on the 

circular rotor lamination due to stator poles 1, 2, 
4, 5, 6, and 8 (N) 

Fj( net lev ^ Total net magnetic levitation force on the 

scalloped rotor lamination due to stator poles 1 , 
2, 4, 5, 6, and 8 (N) 

F t (net) Total net downward magnetic radial force on 

rotor (N) 

Fu net Net downward magnetic radial force on the 
unaligned scallop segment (N) 

F t (u_net) Total net downward magnetic radial force on 

unaligned rotor (N) 

Ft (upper) Magnetic force on circular rotor lamination due 

to stator poles 5, 4, and 6 (N) 

Ft (upper) Magnetic force on scalloped rotor lamination 

due to stator poles 5,4, and 6 (N) 

g General gap length between stator pole and 

rotor (m) 

g a Average stator-rotor gap length at stator pole 1 

or 5 (m) 

g b Gap length between a comer edge of the salient 

stator pole 1 and either of its adjacent rotor 
salient pole comer edges (m) 

g c Gap length between the midpoint of pole face 1 

or 5 and the base of its adjacent scallop cavity 

(m) 

g k Gap length at stator pole k (m) 

go Nominal gap length, 5xl0 -4 m 

i Current (A) 

i k Current in coil k (A) 

L c Circular rotor lamination stack length, 

5.88xl0“ 3 m 

L s Scalloped rotor lamination stack length, 

1.94xl0“ 2 m 

Itw Tooth width, 1.03x10 2 m 

n Number of turns in coil, 215 

P k Permeance at stator pole k (H) 

P k Permeance between stator pole k and the 

circular rotor lamination stack (H) 

Pi Permeance between stator pole k and the 

scalloped rotor lamination stack (H) 
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P T Sum of permeances (H) 

Pj Sum of permeances at the circular laminations 

(H) 

Pj Sum of permeances at the scalloped laminations 

(H) 

R c k Reluctance between stator pole k and circular 
lamination stack (1/H) 

R s k Reluctance between stator pole k and scalloped 

lamination stack (1/H) 

Rj Total reluctance at circular rotor lamination 

stack (1/H) 

Rj Total reluctance at aligned rotor scalloped 

lamination stack (1/H) 

Rj X Total reluctance at scallop cavity through pole 1 

(fig. 10) (1/H) 

Rq Nominal reluctance at the scalloped segment 

(1/H) 

R5 ratio Reluctance ratio at pole 5 


Rfi Parallel reluctance at circular rotor lamination 

stack (1/H) 

R|[ Parallel reluctance at scalloped lamination stack 

(1/H) 

r Outer radius of rotor laminations (m) 

T Torque 

V Flux volume between stator pole and rotor (m 3 ) 

w Magnetic energy (J) 

x Rotor vertical displacement from nominal 

position (m) 

0 Rotor angle, degrees 

§ c k Magnetic flux between stator pole k and circular 

laminations (Wb) 

§1 Magnetic flux between stator pole k and 

scalloped laminations (Wb) 

<1)* Magnetic flux at stator pole k (Wb) 

p 0 Permeability of free space, 47ixl0“ 7 H/m 
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